Entanglement in quantum critical phenomena 
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Quantum phase transitions occur at zero temperature and involve the appearance of long-range 
correlations. These correlations are not due to thermal fluctuations but to the intricate structure 
of a strongly entangled ground state of the system. We present a microscopic computation of the 
scaling properties of the ground-state entanglement in several 1D spin chain models both near and 
at the quantum critical regimes. We quantify entanglement by using the entropy of the ground state 
when the system is traced down to L spins. This entropy is seen to scale logarithmically with L, 
with a coefBcient that corresponds to the central charge associated to the conformal theory that 
describes the universal properties of the quantum phase transition. Thus we show that entangle- 
ment, a key concept of quantum information science, obeys universal scaling laws as dictated by the 
representations of the conformal group and its classification motivated by string theory. This con- 
nection unveils a monotonicity law for ground-state entanglement along the renormalization group 
fiow. We also identify a majorization rule possibly associated to conformal invariance and apply the 
present results to intèrpret the breakdown of density matrix renormalization group techniques near 
a critical point. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.67.Hk 
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The study of entanglement in composite systems is 
one of the major goals of quantum information science 
I], |j, where entangled states are regarded as a valuable 
resource for processing information in novel ways. For in- 
stance, the entanglement between systems A and B in a 
joint pure state \^ab) can be used, together with a classi- 
cal channel, to teleport or send quantum information || . 
From this resource-oriented perspective, the entropy of 
entanglement ab) measures the entanglement con- 
tained in \^ab) (§]■ It is defined as the von Neumann 
entropy of the reduced density matrix pA (equivalently 
Pb), 



E(^ab) = -tr (pa log 2 Pa) 



(D 



and directly determines, among other aspects, how much 
quantum information can be teleported by using \^ab)- 
On the other hand, entanglement is appointed to play 
a central role in the study of strongly correlated quan- 
tum systems j| [| 0, since a highly entangled ground 
state is at the heart of a large variety of collective quan- 
tum phenomena. Milestone examples are the entangled 
ground states used to explain superconductivity and the 
fractional quantum Hall effect, namely the BCS ansatz 
H and the Laughlin ansatz Ground-state entangle- 
ment is, most promisingly, also a key factor to understand 
quantum phase transitions |ï^, O], where it is directly 
responsible for the appearance of long-range correlations. 
Consequently, a gain of insight into phenomena includ- 
ing, among others, Mott insulator-superfluid transitions, 
quantum magnet-paramagnet transitions and phase tran- 
sitions in a Fermi liquid is expected by studying the 
structure of entanglement in the corresponding under- 
lying ground states. 



In the following we analyze the ground-state entangle- 
ment near and at a quantum critical point in a series of 
1D spin-1/2 chain models. In particular, we consider the 
hamiltonians 
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which contain both first-neighbor interactions and an ex- 
ternal magnètic field, and are used to describe a range of 
1D quantum systems [jïo). Notice that Hxy{i — 1) cor- 
responds to the Ising chain, whereas Hxxz(A = 1) de- 
scribes spins with isotropic Heisenberg interaction. Both 
Hamiltonians coincide for 7 = A = 0, where they become 
the XX model. 

Osterloh et al pl| and Osborne and Nielsen O] have 
recently considered the entanglement in the XY spin 
model, Eq. (||), in the neighborhood of a quantum 
phase transition. Their analysis, focused on single- 
spin entropies [ fïï| and on two-spín quantum correlations 
JÏT| , [ï^| , suggestively shows that these one- and two- 
spin entanglement measures are picked either near or at 
the critical point. Here, alternatively, we undertake the 
study of the entanglement between a block of L contigu- 
ous spins and the rest of the chain, when the spin chain 
is in its ground state \^f g ). Thus, the aim in the fol- 
lowing is to compute the entropy of entanglement, Eq. 
(0), for the state \^g) according to bipartite partitions 
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parameterized by L, 



S L = -tr(pL log 2 PL), 



(4) 



where pL — trg |\E , g )(\I r 9 | i s the reduced density matrix 
for Bl, a block of L spins. 

The XY model, Eq. (Q), is an exactly solvable model, 
in that Hxy can be diagonalized by first using a Jordan- 
Wigner transformation into fermionic modes and by sub- 
sequently concatenating a Fourier transformation and a 
Bogoliubov transformation (see for instance |ÏQ| ). The 
calculation of Sl, as sketched below, also uses the fact 
that the ground state \^f g ) of Hxy and the corresponding 
density matrices pl are all gaussian states that can be 
completely characterized by means of certain correlation 
matrix of second moments. 

More specifically, let us introduce two Majorana oper- 
ators, C21 and C21+1, on each site / = 0, • • • , N — 1 of the 
spin chain, 
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Operators c m are hermitian and obey the anti- 
commutation relations {c m ,c n } — 2S mn . Hamiltonian 
Hxy can be rewritten as 



N-l 
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C21+1C21+2- 



(1 - 7) C21C21+3 ] + C21C21+1 ) , 



and can be subsequently diagonalized by canonically 
transforming the operators c m . The expectation value 
of c m when the system is in the ground state, i.e. (c m ) = 
glcml^ g) , vanishes for all m due to the Z2 symmetry 
(of ,af ,af) — » (—af,—af,af) VZ of Hxy- In turn, the 
expectation vàlues 



(6) 



completely characterize \^ g ), for any other expectation 
value can be expressed, by using Wick's theorem, in 
terms of (c m c n ). Matrix B reads [pl| 



B = 
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with real coefhcients gi given, when N — > cxd, by 



1 



91 = TT 
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a cos 4> — 1 — iaj sin < 



2-7T J \acos(f> — 1 — ia7sin< 



(7) 



(8) 



From Eqs. (|7])-(||) we can extract the entropy Sl of Eq. 
(^) as follows. First we compute the correlation matrix of 



the state pl for block £>l, namely S mn + i(Bi,) mn , where 

n n x ••• n L _i" 



b l = 



n 



l-L 
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is constructed by eliminating 2{N — L) contiguous 
columns and rows from B, those corresponding to the 
N — L traced-out spins. Let V G SO(2L) denote an 
orthogonal matrix that brings Bl into a block-diagonal 
form, that is 



B L = VBlV 1 = 



•.L-l , 



1 

-1 



u m > 0. (10) 



Then the set of 2L Majorana operators d m = 
Snto 1 v ninCn, obeying {d m ,d n } = 2S mn , have a block- 
diagonal correlation matrix (d m d n ) = S rnn + i(BL)mn- 
Therefore, the L fermionic operators bi = (d2i+ic?2i+i)/2, 
obeying {ò m ,ò„} = and {b\ n b n } = S rnn , have expecta- 
tion vàlues 



(ò m ) = 0, (b m b n ) = 0, (b\ n b n ) = S„ 



1 + f„ 



fil 



Eq. ( |TT| ) indicates that the above fermionic modes are in 
a product or uncorrelated state, that is 



Pl = Qo 



Ql-i, 



(12) 



where g m denotes the mixed state of mode m. The en- 
tropy oï pl is a sum over the entropy -í?2((l + ^m)/2) of 
each mode [H 2 (x) = —xlog 2 x — (l — x) log 2 (l — x) is the 
binary entropy] and thus reads 
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L-l 
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For arbitrary vàlues of (0,7) in Hxy and in the thcr- 
modynamic limit, N — > 00, one can evaluate Eq. (|^) 
numerically, diagonalize Bl in Eq. ([)]) to obtain v m 
and then evaluate Eq. (|Ï3|). The computational effort 
grows only polynomically with L and produces reliable 
vàlues of Sl for blocks with up to several tens of spins. 
However, further analytical characterization is possible in 
some cases, which speeds the computation significantly. 
For instance, for the XX model with magnètic field, 
7 = and 1/a G [—1,1], one obtains go = (t> a /^ — 2 
and gi — 2sm{l4> a )/l'K for / ^ 0, where <f> a = arccos(l/a), 
from which Sl can be numerically determined. We ob- 
tain, up to L = 100 spins, 



S 



XX 



•log 2 (L) + fci(a), 



(14) 



where fci(a) depends only on a. For the Ising model with 
critical magnètic field, 7 = 1, a = 1, one finds gi = —2/1 
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for odd / and gi 
spins, 



for even l and, for up to L = 100 



k 2 « tt/3. (15) 



Finally, for the Ising model with magnètic field, 7=1, 
and for a close to 1, Kitaev has obtained an analyti- 
cal expression for the entropy of half of an infinite chain 
JÍ5| . To use this result in our setting we need to double 
its value because the entropy resides near the boundary. 
Thus we get 



c<Isinq ~ 1 



1 



1 



(16) 



The XXZ model, Eq. (Q), cannot be analyzed using 
the previous method. Instead we have used the Bethe 
ansatz jïfj to cxactly determine, through a numerical 

procedure, the ground state \^g^) of Hxxz for a chain 
of up to N = 20 spins, from which S'j 20 '' can be computed. 
Recali that in the XXZ model the non-analyticity of the 
ground-state energy characterizing a phase transition oc- 
curs already for a hnite number TV of spins in the chain, 
since it is due to level-crossing. It turns out that, cor- 
respondingly, already for N = 20 spins one can observe 
a distinct, characteristic behavior of S^ 20 ^ depending on 
whcther the vàlues (A, A) in Eq. (||) belong to a critical 
regime. 

The results of the computation of Sl for the spin chains 
(||) and (||) can be summarized as follows. 

Non-crüical regime. For those vàlues (a, 7) or (A, A) 
for which the models are non-critical, the entropy of en- 
tanglcmcnt Sl either vanishes for all L [e.g. when a 
sufhciently strong magnètic field aligns all spins into a 
product, unentangled state] or grows monotonically as a 
function of L until it reaches a saturation value S max . 
For instance, in the infinite Ising chain the saturation 
entropy S max is given by Eq. (|ï^). As shown in Fig. (Q), 
Sl often approaches S max already for a small number L 
of spins. 

Critical regime. Instead, critical ground-states are 
characterized by an entropy Sl that diverges logarith- 
mically with L, 



Sl 



c + c 



log 2 (L) + fc, 



(17) 



with a coefficient given by the holomorphic and antiholo- 
morphic central charges c and c of the conformal field the- 
ory that describes the universal properties of the phase 
transition p7[ , see Fig. (||). This expression was de- 
rived by Holzhey, Larsen and Wilczek [pl for the ge- 
omètric entropy (analogous of Eq. (Q) for a conformal 
field theory), and our calculation confirms it for sev- 
eral critical spin chains. Thus, the critical Ising model 
corresponds to a free fermionic field theory, with cen- 
tral charges c/ = Cf — 1/2, whereas the rest of critical 




Figure 1: Non-critical Sl for the Ising model, Hxy{"/ = 1) in 
Eq. as a function of the size L of the spin block (left axis) 
and parameter a (right axis). The finite correlation length 
governing non-critical phenomena translates into a finite en- 
tanglement length, that is, a finite value of L for which adding 
new spins to a block does not increase its entanglement with 
the rest of the chain. Such an entanglement length (which 
scales as the correlation length) diverges only at the critical 
point a = 1. Scaling arguments also imply that the entropy 
surface is given by log[L/(L|l — a\)]. For any give n a, the 
saturation value for the entropy is given by Eq. (ILa). 



regimes in (||) and (|^) are described by a free bosonic field 
theory, c& = q, = 1 (cf. Eqs. (|l·|) and (|Ï5|)). In particu- 
lar, the marginal deformation < 7 < 1 for a = 1 shows 
scaling for every 7 with universal coefficient c = c = 1/2. 
It is possible to compare the subleading correction be- 
tween two diffcrent vàlues of 7. The behavior we obtain 
is described by 

hm [S L (l = 1) - S L (rf)] = -^log 27 ■ (18) 

The singular behavior at 7 = is the signature of the 
fact that that point belongs to the abrupt change of uni- 
versality class for the XX model. 

The above characterizations motivate a number of ob- 
servations, that we move to discuss. 

Critical and non-critical ground states contain struc- 
turally different forms of quantum correlations. Non- 
critical ground-state entanglement corresponds to a 
weak, semi-local form of entanglement driven by the ap- 
pearance of a length scale, e.g. a mass gap. Indeed, for 
any L, the reduced density matrix çl is supported on a 
small, bounded subspace of the Hilbert space of the L 
spins, and can be obtained by diagonalizing the hamilto- 
nian corresponding to the block Bl and only a few extra 
neighboring spins, as skillfully exploited in White's den- 
sity matrix renormalization group (DMRG) techniqucs 
p9[ . We note here that a bounded rank for pl (relat- 
edly, a saturation value for Sl) is instrumental for the 
success of DMRG schemes, where only of a finite number 
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Figure 2: Numerical calculation of critical entanglement Sl 
for 1D spin chains. Points (+) come from an infinite Ising 
chain with critical magnètic field, Hxy(a = 1,7 = 1) in 
Eq. (g), and corresponds to S L ain3 « log 2 (L)/6 + k 2 . Curve 
(x) comes from an infinite XX spin chain without magnètic 
field, Hxy{cl = 00,7 = 0) in Eq. (g), and corresponds to 
~ log 2 (í/)/3 + fci . Thus the growth in L of entangle- 
ment in the XX model, asymptotically described by a free bo- 
son, is twice that in the Ising model, corresponding to a free 
fermion, Sffj - Sf x M 2(5 , [ s ;™ 3 - S L sinB ). Finally, curve 
(*) comes from a XXX chain of 20 spins without magnètic 
field, Hxxz{A = 1,A = 0) in Eq. (g). These finite-chain 
results, obtained using the Bethe ansatz, combine the loga- 
rithmic behavior (L = 1, ■ ■ ■ , 5) of a free boson field theory 
with a finite-size saturation effect (L = 6, • ■ ■ , 10). We have 
added the lines [log 2 (L) + tt] both for bosons and fermions 
to highlight their remarkable agreement with the numerical 
diagonalization . 



of eigenvectors of pl can be kept. Critical bound-state 
entanglement corresponds, on the contrary, to a stronger 
form of entanglement, one that embraces the system at 
all length scales simultaneously. DMRG techniques have 
reportedly failed to reproduce quantum critical behavior 
and we may, in view of Eq. (|Ï7|), be in a position to 
understand why. Indeed, the divergent character of Sl 
is just one particular manifestation of the fact that the 
number of relevant eigenvectors of pl unboundedly grows 
with L. If, as is the case in DMRG schemes, only a finite 
number of levels can be considered, then a sumciently 
large L will always make the computation of pl impos- 
sible (arguably, even in an approximate sense) by using 
such schemes. This strongly suggests that overcoming 
the above difhculties necessarily requires techniques that 
do not attempt to reproduce the critical behavior of the 
ground state through a local, real space construction. 

Another remarkable, far-reaching fact is that, as men- 
tioned below Eq. (^7|), our results coincides with en- 
tropy computations performed in conformal field the- 
ory. There, a geomètric or fine-grained entropy analo- 
gous to Eq. (||) but for a continuous field theory has 
been considered by sever al authors, including Srednicki 



H|, Callan and Wilczek [||, Holzhey et al [£§ and Fi- 
ola et al |23fl . Thus, starting from non-relativistic spin 
chain models, and by performing a microscopic analysis 
of a relevant quantity in quantum information, we have 
obtained a universal scaling law for entanglement that is 
in full agreement with previous findings in the context 
of, say, black-hole thermodynamics in 1+1 dimensions 

The above connection has a number of implications to 
be exploited. For instance, Srednicki pï[ has obtained 
the behavior of entropy in 2+1 and 3+1 dimensional con- 
formal field theories. For a region 1Z in 2 or 3 spatial 
dimensions, the entropy of TZ is proportional to the size 
o~(lZ) of its boundary, 

S n sa ku{TZ). (19) 

That is, the entropy per unit of boundary area, K, is inde- 
pendent of the size of 1Z. [This is in sharp contrast with 
the same quantity in 1D, where the boundary consists of 
two points and S-ji/a{TZ) diverges logarithmically with 
the length L of 7Z] . Accordingly, Eq. (^9|) also describes 
the critical ground-state entanglement of 2D and 3D spin 
lattices. 

Also the fact that the entropy of entanglement for ÍD 
critical spin chains, Eq. (^7|), matches well-known con- 
formal field theory parameters carries an extra bonus. 
The coefficient in control of the divergent behavior of Sl 
at critical points is the central charge, which is subject 
to Zamolodchikov's c-theorem . The c-theorem states 
that the central charges associated to the ultraviolet and 
infrared end points of renormalization group flows, la- 
beled by Cuv and Cm, obey the inequality Cjjv > G'ir 
for unitary theories. This powerful result establishes an 
irreversible arrow as renormalization group transforma- 
tions are performed. The translation of this idea to the 
quantum information setting is that entanglement de- 
creases along renormalization group flows. An infrared 
theory carries less global entanglement than the ultravi- 
olet theory where it flowed from. The c-theorem seems 
natural as renormalization group transformations inte- 
grate out short distance degrees of freedom, accompanied 
with their quantum correlations. Yet, it is not at all triv- 
ial due to, first, the infinite degrees of freedom (needing 
regularization) existing in a quantum field theory and, 
second, the rescaling step in the renormalization group 
transformation. It is noteworthy, then, that entangle- 
ment decreases both (i) under the local operations and 
classical communication and (ii) along renormalization 
group trajectòries. The former case corresponds to local 
manipulation of an entangled system while the second is 
made out of a block-spin transformation followed by a 
rescaling of the system. Both actions do reduce quantum 
correlations and become irreversible pH p5[ . 

One more remark. From Eqs. (||)-(p|) the complete 
spectrum of Pl can be extracted. The 2 eigenvalues are 
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given by 



L-l 

n 



i + (-ir 



= 0, 1 Vm. (20) 



This allows us to look in more detail to the reshuffling 
of the ground state as more sites are incorporated in the 
block Bl- Every time a new spin is added, the amount of 
local surprise due to quantum correlations with the rest of 
the chain increases, and so does the entropy. But critical 
quantum correlations entangle every single subset of the 
system, and the way they are reordered is far more subtle 
than the relation hinted by entropy arguments. We have 
numerically verified that also a majorization relation ]2^| 
holds for the ground-state reduced density matrices of 
the infmite XY spin chain, namcly 



A 



L+2 



<Xl VL, 



(21) 



where the jump in steps of two is forced by the subtleties 
of the microscopic model. Thus, a critical ground state 
orderly redistributes weights, so as to accommodate for 
the new correlations, according to a detailed, exponcn- 
tially large set of inequalities as contained in Eq. (pi]). In 
this sense, majorization may be a signature -admittedly 
a very refined one- of conformal invariancc. 

But the majorization counterpart in the continuum 
conformal field theory is not yet known. Perhaps, then, 
the correspondence between concepts of quantum infor- 
mation science and conformal field theory, between crit- 
ical ground-state entanglement and geometrical entropy, 
can also be exploited in the reverse direction. 
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